Spin Polarized Ground State for Interacting Electrons in Two Dimensions 
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We study numerically the ground state magnetization for clusters of interacting electrons in two 
dimensions in the regime where the single particle wavefunctions are localized by disorder. It is 
found that the Coulomb interaction leads to a spontaneous ground state magnetization. For a 
constant electronic density, the total spin increases linearly with the number of particles, suggesting 
a ferromagnetic ground state in the thermodynamic limit. The magnetization is suppressed when 
the single particle states become delocalized. 

PACS numbers: 71.10.-w, 75.10.-b, 75.50.Lk 



Ferromagnetic instabilities result from the interplay 
between the electronic Coulomb interaction and the Pauli 
principle. In the Pauli picture, electrons populate the 
non-interacting orbitals of a system, such as a quantum 
dot or a metallic grain, in a sequence of spin up - spin 
down electrons. The resulting minimum spin state min- 
imizes the kinetic energy: it costs energy to flip a spin 
since it must be promoted to a higher energy level. Thus 
the total spin of the system is 5 = when the number 
of electrons N is even and S = 1/2 at odd N. In con- 
trast, the maximum spin allows a maximally antisym- 
metric coordinate wavefunction, thus reducing the effect 
of the Coulomb repulsion (a familiar example of this is 
Hund's rule for atoms). This leads to the Stoner insta- 
bility Q , which gives a spontaneous magnetization when 
the typical interaction exchange energy between two par- 
ticles close to the Fermi level is of the order of the single 
particle level spacing. 

Spontaneous ground state magnetization gives rise to 
interesting effects, which are in the focus of many recent 
studies. In quantum dots, a ground state spin polar- 
ization can explain the absence of an even-odd asym- 
metry in the addition spectra in the Coulomb blockade 
regime Also, the addition of an electron to the 

dot may flip the spin of other electrons already in the 
dot: if the total spins of the ground states of successive 
number of electrons differ by more than 1/2, spin selec- 
tion rules suppress the corresponding conductance peak 
(spin blockade) Spontaneous magnetization effects 
could also explain the presence of kinks in the magnetic 
field dependence of the Coulomb blockade peak positions 
[^,0. The stability of the minimum spin ground state 
in a quantum dot was analyzed for weak interactions in 

. Within perturbation theory, the effective interaction 
strength is enhanced by the presence of disorder, leading 
to a ferromagnetic instability already below the Stoner 
threshold In the diffusive regime, recent studies have 
also considered the effect of mesoscopic wave function 
fluctuations pO] and of off-diagonal interaction matrix 



elements beyond the mean fleld treatment |T|]. The ap- 
pearance of local magnetic moments has been discussed 
also in the strongly correlated limit, at the quantum melt- 
ing of the Wigner crystal |l2[ . 

Although the Stoner instability signals the presence of 
short range magnetic ordering, it is not clear if it will 
also lead to a ferromagnetic ground state in the ther- 
modynamic limit. Actually, the Stoner criterion is ob- 
tained within the mean field Hartree-Fock approximation 
and overestimates the long-range magnetic ordering, pre- 
dicted also in one and two dimensions for the Hubbard 
model at finite temperatures, thus violating the Mermin 
and Wagner theorem j^]. The possibility of a ferromag- 
netic phase in strongly correlated two-dimensional (2D) 
systems was considered in ||lj,|l^ (see also Refs. JTsIp^ ) 
and has recently received experimental support in dilute 
2D electron gases |l^]. At the same time, recent studies 
of fermionic models with random two-body interactions 
show that the ground state polarization is strongly re- 
duced by off-diagonal interaction matrix elements ||ll|,|l8| . 

In this Letter, we investigate numerically the possibil- 
ity of a ferromagnetic ground state in the regime where 
the single particle localization length is smaller than the 
system size. Without interaction this condition is always 
satisfied in two dimensions in the limit of large system 
size !§. 

We study a disordered square lattice with N fermions 
on sites. The Hamiltonian is defined by 



<i,j>,(T iff 
i i^},cr,a' 



(1) 



where C;^ (cia) creates (destroys) an electron at site i 
with spin a, nia = c\^Cia is the corresponding occupation 
number, the hopping term V between nearest neighbors 
characterizes the kinetic energy, random site energies Ci 
are taken from a box distribution over [—W/2, W/2], Ur 
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and U measure the strength of the on-site Hubbard inter- 
action and of the Coulomb interaction, respectively, and 
|i — j| is the inter-particle shortest distance computed on 
a 2D torus (periodic boundary conditions are taken in 
both directions). In the following we choose Uh — U. 

The Hamiltonian (|^) commutes with the total spin 
{[S'^,H] — 0) and its component along an arbitrary 
z-direction {[Sz,H] — 0). Therefore H can be writ- 
ten in a block-diagonal form, with + 1 blocks where 
Sz = -N/2, -N/2 + 1, iV/2, respectively. We consider 
the block with Sz = only, since it is sufficient for analy- 
sis of the ground state magnetization. Indeed, due to spin 
rotational symmetry, the system has a 25* -I- 1 degeneracy 
(where S is the total spin), with Sz = —S, +S; there- 
fore all the eigenvalues of the Hamiltonian (|l|) belong to 
the spectrum of the Sz — subspace. 

The numerical studies of the model ([|) at a finite den- 
sity of interacting particles above a frozen Fermi sea are 
performed in the following way: 

(i) Single particle eigenvalues ea and eigenstates (or- 
bitals) (l>a{i) {ct = 1,...,L^) are obtained via numerical 
diagonalization of the Hamiltonian at Uh ^ U ^ 0. 

(ii) The Hamiltonian (|l|) is written in the basis of non- 
interacting orbitals obtained in (i): 

U E ClldlJ^.'ds.'d,a, (2) 

with — 0Q,(i)cj^, and transition matrix elements 
^2^ = E<^"W'^/3«'^7(i)0*(i), (3) 
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(t>a{i)(t>l3{i)(t>j{i)(t>5ij) 



(4) 



(iii) The Fermi sea is introduced by restricting the sums 
in (H) to orbitals with energies above the Fermi en- 
ergy emf- ct;/5;7:'^ > ^p. We cousidcr a filling fac- 
tor vp = Mp/L/^ = 1/4 (corresponding to 2Mp frozen 
electrons due to spin degeneracy) and a finite density 
p = N/l? of N interacting particles above the Fermi 
level. The frozen Fermi sea approximation is introduced 
for the sake of simplicity, since it allows us to avoid the 
band tail, where the single particle density of states and 
the one body localization length have a strong energy de- 
pendence. The advantages of such an approach have been 
demonstrated in |^ . However, we have also checked that 
the results presented in this Letter are qualitatively sim- 
ilar when Mp = 0. 

(iv) The basis of the Slater determinants, built from 
the single particle orbitals 4'a, is energetically trun- 
cated at high energy orbitals by means of the condition 



Ei=lim^-Mp) < M. Here rui is the orbital index for 
the i-th quasiparticle (to^ > Mp). The truncated Hamil- 
tonian still commutes with the total spin when 5*2 = 0. 
(v) We diagonalize the many-body truncated Hamilto- 
nian. Then the total spin of a given eigenstates is 
found via the application of the operator 5^: S"^]^/) — 
Si{Si + (we take fi = 1). 

We consider iV = 2, 4, 6, 8 particles on a square lattice 
of size L = 8, 11, 14, 16 respectively, at an approximately 
constant density p = N/L^ « 1/32, 0.5 < U /V < 2, 
2 < W/V < 10. In particular, we focus on the localized 
regime W = lOV, where the single particle localization 
length Zi « 4 < L. Data are averaged over a number of 
disorder configurations between 200 and 5000. 

The distribution of the energy differences 5E between 
the ground state energies Eo{S) in spin sectors S = 1 
and iS* = is shown in Fig. |l| for iV = 8 particles in 
the localized regime with W = lOV [5E = Eo{S = 
1) — E[){S — 0)). One can see the mesoscopic Stoner 
mechanism: electron-electron interactions give a sponta- 
neous magnetization {5E < 0), with a probability to have 
a polarized state increasing with the growth of interac- 
tion. The inset of Fig. |l| demonstrates that our results 
are stable when the size Nh of the truncated Hilbert 
space is changed by a factor of five. Even though we can- 
not exclude the existence of very slow iV//-variations this 
check shows that the truncation does not significantly 
affect the ground state polarization. 
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FIG. 1. Normalized distribution of the energy differences 
SE between the ground state energies Eo{S) in the spin sec- 
tors S = 1 &nd S = 0. Here SE = Eo^S = 1) - Eols = 0) 
and < \SE\ > is the absolute value of SE averaged over disor- 
der configurations. Data are shown for N = 8 particles on a 
square lattice of size L — 16, disorder strength W = lOV , and 
interaction strengths U = 0.51^ (circles), U = V (squares), 
and U — 2V (diamonds). Inset shows the distribution at 
U = 2V as a function of the size of the truncated Hilbert 
space: Nh = 932 (triangles up), 2097 (diamonds), and 4354 
(triangles left). 
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The main result of our Letter is shown in Fig. |2[ when 
the number of particle is increased (at a constant elec- 
tronic density) the ground state can be found with high 
probability at larger and larger spin values. In the in- 
set one can see that the ground state average magne- 
tization increases linearly with the number of particles, 
< S >« a{U)N with the slope a{U) growing with 
which determines the strength of the interaction ex- 
change term. The extrapolation of the results presented 
in this figure would give a ferromagnetic ground state in 
the thermodynamic limit. 



jrowth with N. 




FIG. 2. Probability distribution p{S) of the ground state 
spin for W = lOV, U = 2V. Here, iV = 2, L = 8 (cir- 
cles), = 4, L = 11 (squares). A'" = 6, L = 14 (diamonds), 
A'^ = 8, L = 16 (triangles). Inset: ground state average spin 
as a function of the number of particles for U — 0.5V (cir- 
cles) and U — 2V (squares). Straight lines give linear fits, 
with < S{U = 0.5) >^ 0.06iV and < S{U = 2) >« 0.347V. 
Here and in the next figures the error-bars show the size of 
statistical errors. 

The dependence of the average magnetization on the 
disorder strength is shown in Figj^. One can see that 
for sufficiently strong interaction disorder favors ground 
state spin polarization. Indeed, for U/V — 0.5 at strong 
disorder, W/V — 10, the total spin remains less than 
0.5 while for U/V = 2 it becomes five time larger. The 
significant average magnetization < S > appears in the 
localized single particle phase {W/V = 7,10), while in 
the delocalized regime {W/V = 2,4) it remains rather 
weak. 

This is further confirmed in Fig. ^, which shows the 
size dependence of the average spin in the delocalized 
regime W = 2V. With the change of the number of parti- 
cles between TV = 2 and = 8 the average magnetization 
< S > remains constant, with non monotonous fluctua- 
tions around its average value. This is in a sharp contrast 
with magnetization behavior in the localized phase (see 
inset in Fig. 2) where < S > demonstrates a monotonous 




W/V 

FIG. 3. Ground state average spin < S > as a function of 
disorder strength W/V, for U = 0.51/ (circles) and TJ = 2V 
(squares); A'^ = 6 (empty symbols) and N = 8 (filled symbols). 
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FIG. 4. Ground state average spin < S > as a function of 
the number of particles, for W = 2V, U = 0.51^ (circles) and 
U = 2V (squares). The system size changes from L — 8 for 
TV = 2 to L = 16 for AT = 8. 

The ensemble of our results allows us to propose the 
following physical scenario. In the regime when the Fermi 
energy is larger than the electron-electron interaction 
the perturbation theory |^ tells that the correction to 
the spin susceptibility dx induced by interaction is in- 
versely proportional to the conductance of the sample 
g: 6x oc 1/g. With the increase of disorder g drops 
and becomes of the order of one when the single parti- 
cle localization length li becomes comparable with the 
sample size. This indicates that the spin effects become 
more important in the regime of strong disorder, that 
is in agreement with our results (see Fig. 3). In the 
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nonperturbative diffusive regime, non-diagonal interac- 
tion matrix elements start to give quantum fluctuations 
beyond mean-field Stoner approach. These interaction 
fluctuations favor small spin values, since the number of 
off-diagonal scattering events is larger in the lower spin 
sectors of the Hilbert space. This effect can prevent the 
ground state spin from achieving a full polarization [ pT| . 
On the contrary, in the localized regime the off-diagonal 
fluctuations are strongly reduced: due to single particle 
exponential localization. Coulomb repulsion can induce 
electron jumps only inside the localization domains, all 
the other scattering events giving an exponentially small 
contribution. Therefore, a possible scenario is the follow- 
ing: Stoner instability gives, at strong enough interaction 
and/or disorder, spin polarization in domains of the size 
of the single particle localization length, then the cou- 
pling between these domains gives global magnetization. 
The long range nature of the Coulomb interaction seems 
to play a crucial role in this physical picture. Indeed, 
recent quantum Monte Carlo studies of the ground state 
magnetization in the Hubbard model with disorder show 
disappearance of any magnetic order at strong disorder 

H- 

In summary, we have shown that Coulomb repulsion 
can lead to spontaneous ground state magnetization. In 
the regime with localized single particle wavefunctions, 
the total spin increases linearly with the number of par- 
ticles. Even though we cannot exclude that this magnetic 
ordering can become limited at some large finite sizes, our 
results suggest the appearance of a ferromagnetic ground 
state induced by disorder and localization in the thermo- 
dynamic limit. 

We thank the IDRIS in Orsay and the CalMiP in 
Toulouse for access to their supercomputers. 
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